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Abstract: We have studied the kinetics of g-deformed 
bosons and fermions, within a semiclassical approach. 
This investigation is realized by introducing a generalized 
exclusion-inclusion principle, intrinsically connected with the 
quantum g-algebra by means of the creation and annihila¬ 
tion operators matrix elements. In this framework, we have 
derived a non-linear Fokker-Planck equation for g-deformed 
bosons and fermions which can be seen as a time evolution 
equation, appropriate to consider non-equilibrium or near¬ 
equilibrium systems in a semiclassical approximation. The 
steady state of this equation reproduces in a simple mode the 
g-oscillators equilibrium statistics. 


In the limit q —> 1, [a;] —> x and one reproduces the 
familiar boson algebra. 

For g-fermion operators b and b', analogously to the 
boson case, we have the anticommutation relations H,|J 

{6, b} = {6 t , 6’*’} = 0, btf + q~ 1 b^b = q~ N . (5) 

The number operator TV, for this case, can take only the 
values n = 0,1 and the above anticommutations imply 
the following relations 

b | 0 >= | 1 >= 0, | 1 >= 6* | 0 > (6) 


PACS number(s): 05.20.-y, 05.30.-d, 05.40.+j, 05.60.+W 


Recently, there has been a great deal of interest in the 
study of the quantum groups and of their applications to 
g-deformed bosons and fermions. Several models [11JG] 
of g-deformed free bosons and fermions gas |p|,[Ipf 1C|] 
have been applied to Bose gas condensation ]il[ . [l2| . 
anyon fields [ fj~3|| , phonon spectrum in 4 TVe ]l^ |, asym¬ 
metric XXZ Heisenberg chain 0 , nuclear and molecular 
physics |T| |Ii|, conformal field theories Q and quan¬ 
tum Lie superalgebras fH| ]. 

The g-boson and g-fermion particles are defined 
by an appropriate deformation of the commutation- 
anticommutation relations of the creation and annihila¬ 
tion operators, modifying the exchange factor between 
permuted particles. 

For g-boson particles, we have the g-deformed commu¬ 
tations nu 

[a, a ] = [a^, a^] = 0, aa) — qa'a = q~ N ; (1) 

the Hilbert space with basis | n > is constructed with 

these features 

a | 0 >= 0, a t | n >= [n + l] 1 ^ 2 | n + 1 > , 
a | n >= [n] 1 / 2 | n — 1 >, TV | n >= n \ n > , (2) 

implying the following operatorial relations 

a)a = [TV], aa* = [1 + TV] . (3) 

The notation [ ] is defined in terms of the deformation 
parameter q as 


*] = 


r - or 


q-q 1 


( 4 ) 


and, as a consequence, we have 

= [TV], = [1 — TV] . (7) 

In the recent past several authors , assuming 

the relations of Eqs.(3) and (7), have derived the statisti¬ 
cal distributions of g-oscillators (let us recall at this point 
that the statistical distributions can not be derived from 
the thermal averages < a^a > or < tfb > because these 
quantities are the average numbers of particles only in 
the undeformed (g = 1) case). 

In this letter we want to show that the thermody¬ 
namic and statistical properties of g-deformed bosons 
and fermions can be derived in the framework of a ki¬ 
netic approach, recently proposed by us The 

crucial point of our model is the introduction in the tran¬ 
sition probability of an exclusion-inclusion Pauli princi¬ 
ple which can be seen as the semiclassical consequence 
of the partial antisymmetrization-symmetrization of the 
quantum wave function. Assuming that the diffusional 
current must satisfy general conditions and that the 
exclusion-inclusion principle be valid in analogy to the 
pure fermion-boson case, we deduce a generalized non¬ 
linear Fokker-Planck equation appropriate to study the 
equilibrium and also the non-equilibrium conditions of 
g-oscillators. 

Starting point of our approach is the Pauli master 
equation for the mean occupation function n(T, v ) in 
terms of the transition probability ir (t, v — > u) from the 
state v to the state u. In the case of the one dimensional 
velocity space and in the limit of nearest-neighbor inter¬ 
action, with infinitesimal transition from the state v to 
the state v + dv, the master equation is 

dn(t,v) 

——— = n(t, v — dv —> v) + 7T(t, v + dv — > v) 

—7 r(f, v —> v — dv) — 7 r(t, v —> v + dv) . (8) 
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The present formalism can be extended to a general di¬ 
mension velocity space, for this generalization and a more 
exhaustive discussion of the above master equation we 
send to Ref. [ p3[ . 

To introduce, from a semiclassical point of view, the 
quantum behavior of a generic g-oscillator boson or 
fermion system, we postulate a transition probability 
with a generalized exclusion-inclusion principle as [ p3|]24| 

7T (t,v—* u) = r(t,v,v — u) (p[n(t,v)\ip[n(t,u)] , (9) 

where r(t , v, v — u) is the transition rate from the state 
v to the state u, <p( n ) and are real functions that 
inhibit or enhance the transition probability from a site 
to another one. In fact, <p(n) is a function depending on 
the occupational distribution at the initial state v which 
must satisfy the condition y>(0) = 0 (if the initial state 
is empty, the transition probability is equal to zero) and 
ip(n) is a function depending on the arrival state and sat¬ 
isfies the condition = 1 (if the arrival state is empty, 
the transition probability is not modified). If we choose 
<p(n) = n and tp(n) = 1 we find the standard classical 
linear kinetics. 

Explicit expressions of the functions <p(n) and ip(n) can 
be used to simulate semiclassically a quantum exclusion- 
inclusion principle in the transition probability. In Ref. 
fi~i| . [25f we have studied, in this framework, the general¬ 
ized Haldane exclusion statistics p6f and the quantum 
extension of the non-extensive Tsallis entropy p?[ . 

Because of the quantum meaning of the functions tp(n) 
and it must exist an intrinsic relation between these 
functions and the quantum algebra of the creation and 
annihilation operators matrix elements. As tp(n) is pro¬ 
portional to the probability of finding in the state v the 
occupational number n and 4’{ n ) is proportional to the 
probability of introducing an extraparticle into a state 
with occupational number n, we can postulate the fol¬ 
lowing relations 

<p(n) oc|< n — 1 | a n \ n >| 2 , (10) 


ip(n) oc|< n + 1 | | n >| 2 , (11) 


implying an important connection between the semiclas¬ 
sical kinetic approach and the quantum algebra. 

If we expand up to the second order in power of dv 
the r.h.s. of Eq.(8), we obtain the non-linear generalized 
Fokker-Planck equation |f23| , p4| 


dn d 
dt dv 





+D 



<p(n) 


dip(n) 

dn 


dn 

dv 


( 12 ) 


Eq.(12) is a continuity equation and the factor into the 
square bracket is the particle current that can be seen as 
the sum of two terms. The first is the drift current j drift , 
given by 

3 drif t(t,v) = (^J + . (13) 

This quantity is proportional to (p(n)il>(n ), hence propor¬ 
tional to the probability 7r(n) that the particle get away 
from the v site. This is an anomalous current because it 
is a non-linear algebraic function of n. Only in the case 
of the standard Fokker-Planck kinetics (i p(n) = n and 
ip(ji) = 1) this current is proportional to n. 

The second term in Eq.(12) is a diffusion like current, 
j diff , proportional to dn/dv and given by 

j d ift (*» v ) = D A ( n ) ^ n ) , (14) 

with 

., , ,, , dtp(n) dif>(n) . 

A(n) -y(n) g ^ . (15) 

We call j anomalous diffusional current because the 
diffusional coefficient T>(t, v, n) depends explicitly on v 
and n(t , v), while the standard diffusional coefficient does 
not depend on n(t,v). 

Let us now limit ourselves to those physical processes 
described by Eq.(12) where the diffusional current j diff is 
a standard current or Fick current with the diffusional co¬ 
efficient independent of the distribution function n(t,v). 
To realize this condition it is necessary to determine 
the functions and ip{ n ) that satisfy the equation 

A(n) = A. 

In the case of a pure bosonic kinetics ]23|], one has 
ip(n) = n and V’( n ) = 1 + n or ip(n) = y>{n) with n = 
1+n. For a pure fermionic kinetics |23| ] one has tp{n) = n, 
= 1 — n and, in analogy with the bosonic notations, 
the fermionic condition can be written as ip(n) = ip(n), 
where now n = 1 — n. 

The condition V , ( n ) = <^(fi) with n = 1 + an holds both 
for the pure bosonic kinetics (er = 1) and for the pure 
fermionic kinetics (a = —1). 

Now, we impose that the above conditions are con¬ 
served also in the case the kinetics be defined by an ar¬ 
bitrary ip(n) so that we limit ourselves to consider the 
physical processes which satisfy the two following gen¬ 
eral conditions 


A(n) = A , 

ip(n) = <p(l + an) . (16) 

After derivation respect to n of the first condition and 
taking into account the second condition in Eqs.(16), we 
easily obtain the following differential equation for the 
function tp(n) 


where the drift J = J(t, v) and the diffusion D = D(t, v) 
are respectively the first and the second order momenta 
of the transition rate. 
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p( 1 + an)p in) — p{n)p (1 + an) = 0 . (17) 

Introducing the auxiliary function A(n) = p ( n)/p{n ), it 
is easy to verify, by taking into account Eq.(17), that this 
function is equal to a constant a 2 £ IR and our problem 
is reduced to the following Cauchy problem 

p ( n ) — a 2 p{n) = 0 , 

P( 0)=0 , (18) 

¥>( 1 ) = 1 • 

For a real the solution of Eq.(18) is given by 


near-equilibrium g-deformed systems with small (aver¬ 
age) quantum fluctuations. These phenomena are ac¬ 
tually very important in several applications from con¬ 
densed matter to cosmological and high energy physics 


In stationary conditions (t —> oo), the particle current 
in the square bracket vanishes and Eq.(22) becomes a 
homogeneous first-order differential equation, easily in¬ 
tegrated as 


[ 1 - 


(23) 


<P(n) 


sinh an 
sinh a 


while for a immaginary: a = id 


<P(n) 


sin dn 
sini? 


(19) 


( 20 ) 


Eqs. (19) and (20) can be written using the definition of 
the symbol [ ] of Eq.(4) as 


p(n) = [n] , 


( 21 ) 


where q = e a (for q real we have the expression (19) and 
for q complex we have Eq.(20)). 

Thus, selecting ijj(n) = p[\ + an) and imposing diffu- 
sional Fick current, we find the function pin) and ip{n) 
for g-particles (g-boson for a = 1 and g-fermion for a = 
— 1). These functions (p(n) = [n] and ip(n) = [1 + an]) 
correspond to the Eqs. (10) and (11), previously postu¬ 
lated, with the quantum g-algebra of Eqs.(l)-(7). 

Now, if we include the g-functions p{n) and ijx{n) in 
Eq.(13), in the case of Brownian particles (J = ■yv, D = 
"f//3m, /3 = 1/kT and 7 is a dimensional constant), we 
obtain the Fokker-Planck equation for g-oscillators 


d n 

~dt 


d_ 

dv 


jv [n] [1 + an] + 7 


kT a dn 
m sinh a dv 


( 22 ) 


This equation describes the time evolution of g-deformed 
bosons and fermions in the velocity space. 

Great interest is actually devoted to study the statisti¬ 
cal distribution and the dynamical evolution toward equi¬ 
librium of several physical systems. However, because of 
the huge complexity of the kinetic transport theory on 
a quantum field theoretical basis, it results very difficult 
to describe the non-equilibrium behavior of physically in¬ 
teresting systems. 

Because the quantum dynamics is equivalent to classi¬ 
cal dynamics with the inclusion of quantum fluctuations, 
that in average coincide with the Brownian fluctuations, 
we conclude that the non-linear kinetic equation ( 22 ) can 
describe very close the quantum dynamics of a system of 
g-oscillators. Consequently, the above kinetic approach 
can result appropriate to describe non-equilibrium or 


where e = (3{E — p), E = \mv 2 is the kinetic energy 
and p is the chemical potential which can be evaluated 
by fixing the number of particles of the system. 

Eq.(23) defines implicitly the statistical distribution n(v) 
in stationary conditions; it can be explicitly derived as 


n = 


- -tanh 1 

logg 



sinh log q 
— a cosh log q 


(24) 


For q real the quantity n can be written as 


_ 1 , ( e e - aq 17 ^ 

U 2 logg ° g \e f - aq° ) 

while for q = e 1 ^ complex as 


(25) 


n = — tan 
d 


-1 


sind 


e e — a cos d 


(26) 


We like to stress that to consider complex values of q can 
be relevant in reproducing the exact interaction between 
particles in several physical systems. In Ref. JTi| it is 
analyzed the g-deformed pairing vibration in nuclei and 
is shown that the real part of q simulates an attractive 
residual interaction between the nucleons and the imag¬ 
inary part of q decreases the binding energy of the pair 
of nucleons. 

If we use as variable the single-particle energy e, it is 
easy to verify that the second order density fluctuation 
(An ) 2 =< n 2 > — < n > 2 can be expressed as || 2 d| 


(An ) 2 



<fi(n) 

ip(n) 



(27) 


or more explicitly for g-deformed bosons and fermions 


—-—sinh a r i n i 

(An ) 2 =- n [1 + an) . (28) 

a 

The thermodynamic relations for g oscillators can be 
derived starting from the explicit expression of the <p(n) 
and ifi(n) functions (see Ref. [|4j for details). The entropy 
density can be expressed as S(n) = flog[ift(n)/<p(n)]dn 
or more explicitly as 


s<»)= / iog!i±f!U 
J N 


(29) 
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In conclusion, we have studied the kinetics of the q- 
oscillators defined in Eqs.(l)-(7) in the framework of a 
semiclassical non-linear approach. The two families of 
g-bosons and g-fermions can be interpreted as Brow¬ 
nian particles (J oc v, D = const). They obey to 
an exclusion-inclusion principle defined in terms of the 
transition probability by the functions ip(n) = [n] and 
ip(n) = [1 + an] which are implicitly related with the 
creation-annihilation operator matrix (see of Eqs.(10) 
and (11))- This relation fixes an important connection 
between quantum algebra and the semiclassical kinetics. 
In this framework, we have derived a generalized Fokkcr 
Planck equation appropriate to study near-equilibrium 
and non-equilibrium (/-oscillator systems. This equation 
is easily integrated in stationary condition and repro¬ 
duces the statistical distributions for (/-deformed bosons 
and fermions both for real and complex q values. 

We would like to thank M. Rego-Monteiro for stimulating 
discussions and for the critical reading of the manuscript. 
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